Exact Landau Levels in Two-Dimensional Electron Systems with Rashba and 
Dresselhaus Spin-Orbit Interactions in a Perpendicular Magnetic Field 
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We study a two-dimensional electron system in the presence of both Rashba and Dresselhaus 
spin-orbit interactions in a perpendicular magnetic field. Defining two suitable boson operators and 
using the unitary transformations we are able to obtain the exact Landau levels in the range of 
all the parameters. When the strengths of the Rashba and Dresselhaus spin-orbit interactions are 
equal, a new analytical solution for the vanishing Zeeman energy is found, where the orbital and 
spin wave functions of electron are separated. It is also shown that in this case the Zeeman and 
spin-orbit splittings are independent of the Landau level index n. Due to the Zeeman energy, new 
crossing between the eigenstates \n, k, s = 1, a > and |n -I- 1, k, s' = — 1, ct' > is produced at certain 
magnetic field for larger Rashba spin-orbit coupling. This degeneracy leads to a resonant spin Hall 
conductance if it happens at the Fermi level. 

PACS numbers: 72.20.My, 73.63.Hs, 75.47.-m 



Recently spintronics, operating the spin of charge car- 
riers in semiconductors, has attracted much attention in 
both theory and experiment due to its potential applica- 
tion for information technology [1-7]. One crucial step in 
spintronics is to produce and efficiently control the spin 
current. The Rashba spin-orbit interaction [8], which 
can be adjusted up to 50% by the application of gate 
voltages, is used to arrive this goal [9]. The Rashba cou- 
pling stems from the lack of structure inversion symme- 
try, and can be strong in semiconductors (InGaAs, Al- 
GaAs) [10, 11]. Besides the Rashba coupling there also 
exists a Dresselhaus spin-orbit interaction in semiconduc- 
tors like InSb/InAlSb [12], which originates from the lack 
of bulk inversion symmetry. 

The coexistence of the Rashba and Dresselhaus spin- 
orbit interactions gives rise to some novel spin-dependent 
phenomena in semiconductors. Recently the spin Hall 
effect in two-dimensional electron system in the absence 
of a magnetic field has been studied [13,14]. When the 
Rashba and Dresselhaus spin-orbit interactions have the 
same strength, the spin state of the wave functions is in- 
dependent of the wave vector. The phenomenon was pro- 
posed to design a nonballistic Spin-Ficld-Effcct Transis- 
tor [15] . The experimental observations of a spin-galvanic 
effect and spin-orbit coupling weak localization effects 
have been used to understand the interplay between dif- 
ferent spin-orbit types and to obtain the ratio of the 
Rashba and Dresselhaus spin-orbit coefficients [16,17]. 

In the presence of a perpendicular magnetic field, the 
Rashba and Dresselhaus spin-orbit interactions couple all 
states in each Landau level, and an analytical solution to 
the full problem has been attempted, but without success 
[18-23] . As a result, the spin Hall conductance in this sys- 
tem does not get understood. In Refs. [20,21], Shen et al. 
investigated transport properties of the electron system 
with a pure Rashba or Dresselhaus coupling. A resonant 
spin Hall conductance is found for the Rashba type due 



to the level crossing induced by the Zeeman energy at the 
Fermi level. It is very interesting to know whether the 
resonant spin Hall conductance still exists in the pres- 
ence of both spin-orbit couplings. When the strengths of 
both spin-orbit interactions are equal, an exact solution 
is obtained in the absence of the Zeeman energy [22-25] . 
However, we shall see below that the energy spectrum 
is incomplete. The approximate energy spectrum is ob- 
tained in Refs. [21-23] by using the perturbation the- 
ory. In the present work, we solve analytically this open 
problem and construct the eigenvalues and eigenstates 
for such a system so that the spin Hall effect could be 
investigated in the future. 

The Hamiltonian for a single electron with spin-^ in a 
plane under a perpendicular magnetic field is given by 
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where 11 = p + |A, ai{i = x, y, z) are the Pauli matrices 
for electron spin, gs is the Lande g- factor, /zb is the Bohr 
magneton. Here we choose the Landau guage A = ySx. 
The spin-orbit interaction has the form 
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where a and /3 represent the Rashba and the Dresselhaus 
spin-orbit couplings, respectively. Note that Px = k is a 
good quantum number due to [px, H] = 0. 

In order to diagonalize the Hamiltonian (1), we intro- 
duce the bosonic annihilation operator bka = [z/ + 

-^(}i-\-iPy^^ \/\ajiaD\'u.a\ and the corresponding creation 
operator fej,^ = (fefcff)^, with the cyclotron radius = 

u„ = f7[l - isgniaRau)], 
Here we point out that 
in the operator hha plays an 
important role in solving the Hamiltonian (1). Differ- 
ent the pure Rashba or Dresselhaus coupling case, the 
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and the orbital index cr = ±1. 
the constant \/\aRaD\ua 
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orbital space of electron is divided into two infinitely di- 
mensional subspaces described by the occupied number 
representations T„ associated with bi~„ and b],^ due to 
this constant term (also sec Rcf. [23]). When ud = aji 
and g = 0, the total Hamiltonian (1) can be diagonalized 
in by a unitary transformation [22,23]. It is surpris- 
ing that in the special case the eigenstates obtained in 



Ref. [22-25] are nothing but the basis vectors of two 
subspaces. This means that the Landau levels decouple. 
However, in these works, another solution where all the 
Landau levels couple each other is omitted. In terms of 
the bosonic operators bka and bj,^, the Hamiltonian (1) 
can be rewritten as H = H-i Hi, and 



= ^la^ka + (VWoaRlu^ - iV2aRa+ + V2aDCT-)bl^ + {^/\a^a^\ul + iV^ana- + V2aDO+)bka - -^9(^z + A, (3) 

where the cyclotron frequency lo = f^, a± = \{(Jx ± i<yy), A = 2\aDaR\ + 5, and g = with me the free electron 
mass. 

We can see from Eq. (3) that a state |to, t> in the basis of r^r is coupled to |m + 1, J.> due to the Dresselhaus 
coupling, which is further coupled to \m + 2, t> due to the Rashba coupling. Therefore a Landau level is coupled to an 
infinite number of other Landau levels. This means that in general the eigenstates for the Hamiltonian (3) in T„ are 
convergent infinite series in terms of the parameters ao, clr and g. After a tedious but straightforward calculation, 
we find that the two-component eigenstate for the nth Landau level with quantum numbers s and a has the form 
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where the spin index s = ±1, A^s is a parameter to be determined below by requiring and to be nonzero, 
4)mka is the eigenstate of the mth Landau level in r^-, i.e. bl^(p„ika = Vm + Icpm+ika, bkafpmka = Vrn4>m-ika and 

< (f>m'ka'\4'mkcr >= Smm'Saa'- When m +00, a"^^ = = 0. 

(i) When |a£)| = \aR\ = a, the eigenvalue for the nth Landau level with s and a is given by 

Ens = fujcns = huj{n + ^ + 2a^ + ^sy/g'^ + 640^). (5) 

We note that Ens=-i ioY g = is nothing but the exact eigenvalue obtained previously [22-25]. Therefore, the 
exact energy spectrum is incomplete due to the lack of the eigenvalue for s = 1 in Eq. (5). Obviously, the energy 
splitting induced by the Zeeman energy and spin-orbit coupling is independent of the Landau level index n, i.e. 
AE = Ens=i - Ens=-i = hw^/g^+6Aa\ 

Because the eigenvalues are independent of the parameter A„s, for simplicity, we choose A„s = 0. For the corre- 
sponding eigenstate, we have the recursion relation 

[cns - m + 2a2 - i(l - g)]aZ = [cns -m + 2a^-^{l+ g)](iZ, 



V^au. [aZ-u + PZ-is] - i^ns - m - 2a2 - - (1 - g)]^- + ^ci^f^Zs + V^^aul[aZ+is + PZ^u] = 0. (6) 

Equations (5) and (6) are derived by demanding the coefficient determinant of and PZ to be zero. The same 
approach is employed in deriving Eqs. (7)-(ll) below. When g = 0, from Eq. (6), we have the solution a^s = —Pn's 
and = PZ = 0(m ^ n) for the nth Landau level with spin index s = — 1 in To-, which recovers the exact 
eigenstate [23]. We would like to stress that new solution obtained here describes the electron with spin index s = 1 
in T„. Note that when g — G, there is no spin polarization in each Landau level, i.e. < n,k,s,a\^az\n,k, s,a >= 0. 
So it is not expected that there exists an out of plane spin current. In this case, the orbital and spin wave functions 
of electron are separated, and the spin part -^(1 sT„)'^ is independent of the wave vector fc, which is similar to that 
observed in the absence of a magnetic field [15]. 
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(ii) When |a£)| ^ \afi\, the reduced eigenvalue is given by 
|az)afl|[2|a£,afl|(l - J - ^A^^] 1 



where s = 1 for n = and the parameter A„s is determined by the following equation 
V\0'DaR\{{\aD\ - \aR\)[2{\aD\ + - |ar.|A^g) - gAns\{ens -n- 
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(7) 



+ aD + aji) 



+(1 + A„,)(|aH| - |ar,|A„s)(a|j -a%+ -g)} = 0. 



(8) 



If jaul = or |aij| = 0, Eq. (8) is automatically satisfied and the eigenvalue (7) becomes that in the presence of a 
pure Rashba or Dresselhaus spin-orbit interaction [8, 12]. When g = 2\a%^ — ajj\, we have the analytical expression 
for the eigenvalue 



ens=n + 2\aDaR\ + isy'[l - 2sgn{a\ - a|,)(|aD| + |ai?|)2]2 + 8n(|aD| - |afl|)2. 



(9) 



Here, we have used A„s = sgn{aj^ — a|j). Note that when aDiaR) aRiflD), the energy spectrum described by Eq. 
(9) is identical to that in Eq. (5) with qd = cir and g = 0. 
The eigenstate for the nth Landau level with s and a is determined by the recursion relation 

[B„s(e„s - n - A + 1 + C„s) + \aDaR\AnsVns]K-U = \/nulBns{Bns - Cns)0^l-, 
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where 



A — 



B„ 



(1 - A„s)(|a_R| - |aj)|A„s) 
\aR\ + \aD\Al, 

\aDaR\{l - Al,) 



\aR\ + |a£)|A2^ 



^g(|afl| - laplAp + 4|az3aK|(|aD| + |afl|)A, 
|aii| + |aD|A2^ 



2{\aR\ + \aD\){\aR\~\aD\Al^)-gA 



\aR\ + |aD|A2^ 
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We note that when an = or a^j = 0, the corresponding 
eigenstates presented in Refs. [8,12,20-25] can be easily 
obtained from Eq. (10). 

Fig. 1 shows several low-lying Landau levels as func- 
tions of \aR\ for \(3\ = 0.2|a|, 0.6|a|, and 1.2|a| and g = 
and 0.1. When jaflj is small, the energy levels Ens=i 
and En+is=-i are approximately degenerate for g = 
[see Fig. 1(a), (c) and (e)]. If g ^ 0, the approximate 
degeneracy is lefted in Fig. 1(b), (d), and (f). However, 
the Landau levels shift in an opposite way for ja/f j > jacj 
and jaflj < jacj. When ja/jj > jauj new crossing between 
|n, fc, s = 1, C7 > and \n + 1, k, s' = —1, a' > pointed by 
arrows is produced in Fig. 1(b) and (d) while there is no 
new mixing for jaijj < jauj in Fig. 1(f). This new degen- 
eracy induced by the Zeeman energy leads a resonant spin 
Hall conductance if it happens at the Fermi level, similar 
to the case in the presence of the pure Rashba spin-orbit 
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interaction [20,21]. It is expected that the resonant spin 
Hall conductance can be applied in device design. The 
degenerate point at the Fermi level is easily controlled 
by the gate voltage, which changes the Rashba coupling 
strength. We note that with increasing the Dresselhaus 
coupling, the resonance is moved to the lower magnetic 
fields. 

In summary, we have analytically derived the complete 
energy spectrum in 2D electron systems in the presence 
of both Rashba and Dresselhaus spin-orbit interactions 
under a perpendicular magnetic field. The eigenvalue 
for the nth Landau level with s and a is determined by 
the parameter A„s, which satisfies the highly nonlinear 
equation (8) for \a]j\ ^ \aj^\. For |a£)| — |afl|, we obtain 
a new analytical solution when the Zeeman energy van- 
ishes. Each Landau level possesses two-fold degeneracy 
because it is independent of the orbital index a. We also 
note that the eigenvalues depend only on the absolute val- 
ues of the Rashba and Dresselhaus spin-orbit coefficients 
rather than their signs. The corresponding eigenstates 
are convergent infinite series except those at \aii\ = |a£)|, 
g — and s — ~l. Once the energy spectrum is known, 
we can use it to calculate physical properties of the sys- 
tem, such as charge and spin Hall conductance. The work 
in that direction is in progress. 
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the University of Houston and by the Robert A. Welch 
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FIG. 1: (Color online) Landau levels in unit of huj as a func- 
tion of the dimensionless parameter \aR\. The parameters 
|a_D| and g are shown in the panels. Solid, s = 1; dot, s — —1. 
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